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ABSTRACT 

Blazars radiate from relativistic jets launched by a supermassive black hole along our line of sight; the 
subclass of FSRQs exhibits broad emission lines, a telltale sign of a gas-rich environment and high accretion 
rate, contrary to the other subclass of the BL Lacertae objects. We show that this dichotomy of the sources in 
physical properties is enhanced in their flaring activity. The BL Lac flares yielded spectral evidence of being 
driven by further acceleration of highly relativistic electrons in the jet. Here we discuss spectral fits of multi-/l 
data concerning strong flares of the two flat spectrum radio quasars 3C 454.3 and 3C 279 recently detected in 
y rays by the AGILE and Fermi satellites. We find that optimal spectral fits are provided by external Compton 
radiation enhanced by increasing production of thermal seed photons by growing accretion. We find such flares 
to trace patterns on the jet power - electron energy plane that diverge from those followed by flaring BL Lacs, 
and discuss why these occur. 

Subject headings: galaxies: active — radiation mechanisms: non-thermal — galaxies: BL Lacertae object — 
quasars: general — quasars: individual (3C 454.3, 3C 279) 



1. INTRODUCTION 

Blazars radiate from narrow, relativistic jets that blaze on 
us when closely aligned with our line of sight. The jets 
are launched by a central supermassive black hole (SMBH, 
masses M # ~ 1O 8 ~ 9 M0) with bulk Lorentz factors T up to 
several tens, and yield beamed non-thermal radiations with 
observed fluxes and frequencies enhanced by aberration and 



Dopp ler effects of Special Relativity ( jBegelman et al.|1984 
Koni gl|1986[|Urry & Padovani|1995[ ). 



The subclass of the BL Lac objects show few signs of radi- 
ation from an accretion disk, and weak or intermittent if any 
emission lines; such conditions indicate little surrounding gas, 
and low current accretion rates on the order of m < 10~ 2 in 



Eddington units (e.g., Ghisellini et al. 2009| ). Their spectra 
are well represented as a smooth spectral energy distribution 
(SED) S v = vF v , featuring two peaks: one in the range of ra- 
dio - soft X-ray frequencies, due to synchrotron emission by 
highly relativistic electrons inside the jet; and another at hard 
X-ray or y-ray energies, due to inverse Compton upscattering 
by the same electrons of the seed photons provided by the syn- 
chrotron emission itself (synchrotron self-Compton radiation, 
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SSC; Jones et al.|1974l |Marscher & Gear|1985HMaraschi et| 
|al.|1992p . 

The other blazar subclass is constituted by the Flat Spec- 
trum Radio Quasars (FSRQs); besides the non-thermal peaks, 
they also feature the strong broad emissi on lines and h igh Big 
Blue Bump typical of quasars (e.g., |Peterson||1997[ ), which 
yield evidence of plenty surrounding gas associated to high 
current accretion rates with m ~ 1. Such an accretion pro- 
vides seed photons from outside the jets, that drive an addi- 
tional contribution to the non-th ermal radiations by the exter- 
nal Compton scatte ring (EC; see |Dermer & Sc hlickeise rfl993| 
Sikora et al.||1994| l. These EC con tribution s ofte n dominate 
the y-ray outputs of FSRQs ( |Maraschi et al.|2001| i. 

Both source types exhibit intermittent episodes of major 
brightening on timescales of days or even shorter, with sub- 
stantial flux surges ("flares") particularly at y-ray energies, as 
recently observed for 3C 454.3 by AGIL E-GRTD (|Pacciani et| 
|al.|2010| > and for 3C 279 by Fenni-LKT ( |Abdo et al.|201CT . 

Here we examine how the flares of the two subclasses differ, 
and relate their differences to specific physical processes in 
the sources. Flares in BL Lacs have already been studied in 



Paggi et al. ( 2009a i and related to current acceleration of the 
emitting electrons; so in this paper we focus on flaring activity 
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of FSRQs. 

2. RADIATION: FSRQS VERSUS BL LACS 

The SEDs of all blazars include a basic SSC process. The 
synchrotron component is produced since the jets contain a 
magnetic field B, highly relativistic electrons with random 
Lorentz factors y up to 10 6 - 10 8 , in addition to non-relativistic 
protons with a common density n within a comoving size 
R; all flow towar d the observer with bulk Lorentz factors 
T ~ 10 - 20 (see |Bottcher|[2007) also |Celotti & Ghiseilini| 
2008). Interacting with seed photons emitted inside or outside 
the jet, the same electrons also radiate by inverse Compton in 
the SSC or EC process, respectively. 

2.1. Spectra from SSC 

We will mainly adopt for the electron energy distributions 
curved, log-parabolic shapes; adding to their value as fitting 
tools for the associated spectral energy distributions, they are 
theoretically based as discussed in detail in Sect. H] 

In fact, our electron distributions will have the form 



(1) 



Here jq is the electron injection Lorentz factor and s is a con- 
stant contribution to the slope. On a log-log plot, such dis- 



tributions appear as parabolic humps with curvature r (Mas- 
|saro et al. 2004); the parameters may vary during flares. 
The energetic content of such an electron population is best 
expressed in terms of the rms adimensional energy y p = 

y 2 N{y)dy/ f N{y)dy) * ; this turns out to also provide the 

peak position for the function y 2 N(y), that will enable us to 
compare BL Lacs and FSRQs energetics on a similar footing. 

On the other hand, the SSC spectra radiated b y such dis- 
tributions are them selves curved (Massaro et al. 2004~{ [Tra^l 
|macere eTa l. 2007 and references therein), and their SEDs 
are effectively represented by the log-parabolic shapes 



S(y) = S, 



(2) 



Here the spectral curvature is marked by b, and the peak 
frequency v p and flux S p ; the two SED peaks closely yield 
the total flux in each component. For the synchrotron com- 
ponent one has a peak frequency v$ oc By 2 , a peak flux 

Ss oc B 2 y 2 nR 3 , and the spectral curvature b s w r/5. For 
the IC component one has a peak frequency vq oc By p , a peak 
flux Sc °c B 2 y p n 2 R 4 , and the spectral curvature be ~ r/10, in 
the Thomson regime. 

From an empirical standpoint, such curved spectra have 
been recognized to provide very go od spectral fits with uni- 
form residuals, from the X-ray (e.g., [Landau et al.|1986{|Mas-| 
|saro et al.|2004"l|Tramacere et al.|200?T Massaro et al.|2008| l 
to the sub-mm band ( |Gonzalez-Nuevo et al.|2010| i. 

2.2. Spectra from EC 

In the FSRQs in particular, the higher energy component 
usually features substantially larger fluxes with respect to the 
synchrotron component, at variance with the yield of the SSC 
process in conditions of photon energy density not exceeding 
the magnetic one.. This is often interpreted in terms of an ad- 
ditional EC component originated from seed photons emitted 



outside the jet, specifically by the broad line region (BLR) or 
by the inner accretion disk (|Dermer et al.|199 2 ; Sikora et al. 



1 994] |Bottcher & Dermer|2002) .' 

The disk emission is widely modeled as a sum of annular 
concentric surfaces radiating locally as a black body ; their 
temperatur es follow along the radius % the profile (see jFrank 
|et al |2002| > 



T(r) = T, 



i - 



vl/4 



1/4 



(3) 



where , M is the accretion rate and 

R, « GM m lc 2 is the radius of the last stable orbit. The in- 
tegrated emission has a peak flux S 4 oc M oc Tj 4 at a peak 
frequency Vd °c M 1/4 oc T,. The BLR instead comprises a 
number of "clouds" intercepting some 10% of the radiation 
from the disk at distances ~ IP 17 " 18 cm, with ionization pa- 
rameter values ~ 10 _1 (see Peterson 1997] >. 



These narrow spectra of external photons interact with the 
electrons in the jet to yield EC radiation, with a SED again fol- 
lowing shapes as represented by Eq. [2] For the photons com- 
ing from the disk to reach the emitting region at a height scale 
D, the EC radiation has a peak flux Sj oc My 2 p R} nD 2 T l , 
• 1 9 1 

a peak frequency oc M* y p T~ , and the spectral curva- 
ture bd ~ r close to the electron distribution curvature. For 
photons coming from the BLR we have a peak flux S<g oc 
M y 2 R 3 n T, a peak frequency vg oc Ms y 2 T, while the spec- 
tral curvature is again b<g * r. 

We recall that the observed frequencies are enhanced by the 
beaming factor 5, while the fluxes are boosted by 6 A , w ith 5 ~ 
2rfor small viewing angles ~ 1/r (Begelma n et al.|1984) >. As 
mentioned above, the EC flux is further multiplied by factors 
r 1 -T 2 depending on the scattering geometry for the disk and 
the BLR photons, respectivel y ( Dermer & Schlickeiser|2002 
IGhisellini & Tavecchio|2009) . 

Examples of these emissions and contributions to the SEDs 
are numerically computed and illustrated in Fig. [T] for both 
the low and high states of a typical FSRQ with M # x 10 8 Mq, 
m * 1, r * 10, R * 10 17 cm, B * 1 G and n * 10 2 enr 3 , and 
10 2 . These will be used in Sect. 3.1 as a guide to focus 
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the main flare drivers, preliminary to our detailed fits in Fig. 

El 

3. FLARES OF 3C 454.3 AND 3C 279 

The dominant process driving the blazar flares can be tested 
from the observed spectral changes compared with the tem- 
plates provided by Fig. [T] 

3.1. Spectral variations 

In several flaring BL Lacs the SEDs have been observed to 
inc reas e both their peak heights and positions (see Paggi et 
al. |2009a) >; in particular, the synchrotron peak considerably 
shifting up in frequency constitutes clear evidence of electron 
acceleration. 

On the other hand, in several FSRQs the dominant and 
growing second spectral peak, along with the nearly stable 
synchrotron frequency, favor the EC process driven by in- 
creasing seed photon flux. This is shown in Fig. JTj-f, in com- 
parison with the alternatives presented in Fig. nfa, b, c, d, 
e. Correspondingly, our fits to the specific SEDs of 3C 454.3 
and 3C 279 are presented in Fig. [2] with model parameters 
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Figure 1. T emplates of the SEDs of a typical FSRQ. The pre-flare condition (blue lines) comprises the five spectral components with their parameters discussed 
in Sect. |2.2| synchrotron emission (full lines), SSC radiation (dashed lines), accretion disk emission (dotted lines), disk EC radiation (dot-dashed lines), and BLR 
EC radiation (double-dot-dashed lines). We also represent flare conditions (red lines), driven by different processes corresponding to an increase by a factor 2 in 
the key parameter: electron energy y p (panel a), magnetic field B (panel b), comoving radius R (panel c), height scale D (panel d), bulk Lorentz factor T (panel 
e), and accretion rate m (panel f). In Table^we give the detailed parameter values for the case (f), that is shown to be most relevant in Fig. [2] 



given in Table [T] Notice t he larger particle d ensities found in 
FSRQs than inBL Lacs (Pag gi et al.|2009b) , consistent with 
different matter contents in the surroundings of the two source 
kinds. 

We stress that for both sources the fits require considerably 
larger values of m in the flares, with limited variations of the 
other parameters; for 3C 279, however, the low X rays in the 
high state imply also a larger distance D of the emitting re- 
gion. Source transparency to pair production process via y-y 
interaction is assured by the optical depths r yr = 0.03 and 
T yy = 0.24 for hv < 1 GeV, in the high state of 3C 454.3 and 
3C 279, respectively. 

So the large flares of the FSRQs 3C 454.3 and 3C 279 ap- 
pear to be dominated by thermal seed photons increasing their 
fluxes much more than the frequencies, as e xpect ed from the 
increasing accretion rates discussed in Sect. |2.1| To wit, the 
EC flux grows strongly in y rays while in the IR-opfical bands 
little or no increase occurs for the synchrotron peak frequency 
(see Fig [2]). 

On the other hand, the modest rise of the observed X-ray 
fluxes during flares constrains any increase of n; so increasing 
seed photons apparently dominate the y-ray surges. 

3.2. Jet Power 

To further probe the physical processes driving a flare, an 
important quantity is provided by the overall energetics. The 
total jet power Lj is comprised of the radiative component 



L r = nD 2 L F/4F 2 , and of the kinetic powers in electrons and 
protons L e + L p = 4nR 2 c 3 n V 1 [ m e (y) + m p ) /3, with a sub- 
dominant magneti c component ( |Celotti &~G hisellini 2008| 
Paggi et al. 12009b) . 

Here F Is the flux observed at the luminosity distanc^] 
Dl(z), and m e and m p are the electron and proton masses. 
One "cold" proton per electron is assumed, so for the FSRQs 
with average electron energies (y) < 10 3 the proton compo- 
nent to the energetics is considerable and hardly variable; this 
steadily feeds the energetics of the radio lobes (ICelotfi et al.l 
|2D0T) . 

3.3. Flare evolutions 

Our main aim here is constituted by flares of bright blazars; 
their historical conditions are effectively represented in the L T 
- y p plane (see Fig|3), where they appea r to be strung on av- 
erage along a "bright blazar str i p" (see |Cavaliere & D'Elia 
[2002) [B5ttcher & Dermer|2002| |GhTsellini et al.|2002) . The 



BL Lac sources reside in the right-hand region, around the 
line where the crossing time t cr - R/c matches the radiative 
cooling time t c oc 1 /yw related to the energy density w in the 
magnetic or radiation field, to yield Lj oc On the other 
hand, the FSRQs lie in the left-hand region, around the line 

1 In the following, we adopt the standard flat cosm ology with £2a = 0.74 
and H = 72 km s" 1 Mpc~' (seelDunkley et al.|2009 . 
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1E-8-, 




1E-8-, 




Figure 2. Our fits to the SEDs of two specific FSRQs. (Top panel), low (blu e lines) and h igh (red lines) state of the FSRQ 3C 454.3 as recently reported by 
|Pacciani et"aT]j2010^ ; (bottom panel), low and high states of the FSRQ 3C 279 jAbdo et al.|2010) . Model parameters are given in Table[T] 



where the cooling time is matched by the acceleration time 
t n cc y/E related to the effective electric field E (see Cavaliere 
& D'Eha|200"2l|Ghiselliiii et al.|2009| , to yield L T cc y~ L . The 
lower/left corner of the plane is populat ed by wea ker sources 
out of our present scope (see Pad ovani et al.|2003) . 

We see from Fig. [3] that on this plane major flares in our 
database move in directions different from the slope of the 
bright blazar strip, into regions where the sources have higher 
luminosities and total powers and faster cooling. In particular, 
the FSRQs move vertically, whereas luminous BL Lacs (e.g., 
S50716+714 an d Mrk421) move in a substantially slanted di- 
rection ( Paggi et al.||2009b] >. We interpret these patterns in 
terms of source structure changing as follows. 

The flares of gas-poor BL Lacs are mainly driven by in- 
creasing y causing stronger synchrotron radiation. Cooling 
is faster, but still can be balanced by shorter crossing times 
through a smaller source size (like is the case for jets with 
inner structure, see Tavecchio & Ghisellini 2008, Giannios et 
al. 2009), to yield y ~ l/R and L r cc N/R 2 in terms of the 
electron number N cc nR^ . Eventually, the source power may 
exceed the yield from the feeble current accretion; so it has 
to live on the alternative reservoir constituted by the the ro- 



tational energy of the SMBH, that was accumulated during 
previous massive accretion episodes, and is now tapped via 
the Gene ral Rel ativistic mechanism proposed by Blanford & 
This implies L < 8 x 10 45 M./10 9 M© erg/s; 



1977 



Znajek in 

weaker and more massive sources like Mrk501 are not af- 



fected by such a constraint (see also discussion by Paggi et 
|al.|2009b| l. 

On the other hand, the flares of the luminous FSRQs are 
not so constrained in view of their high current accretion; the 
sources increase their luminosity L r cc y 2 NS c /, and in enter- 
ing the fast cooling region are assisted by acceleration from 
E increasing with luminosity to yield y cc R y/E/L r ~ const 
(see [Cavaliere & Morrison 198 ). Note that L T > L r holds af- 
ter the discussion in Sect. 3.2 implying on the Lj - y p plane 
relatively short excursions even during substantial radiative 
flares; we have extended the relative arrows in Fig. [3]to high- 
light their directions, the relevant information in the present 
context. 

Thus in our synoptic view of Fig. [3] a relation is apparent 
between the source location and the flaring patterns of blazars 
on the Lj - y plane. 
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4. DISCUSSION AND CONCLUSIONS 

Before concluding, we discuss why such a relation is not 
biased by our choice of log-parabolic shapes for the elec- 
tron energy distributions and related SEDs given in Eqs. [T] 
and [2] We note that cooling (particularly in conditions of 
bright EC radiation) will erode them at high energies from 
their basic shapes. Actually, numerical simulations includ- 
ing synchrotron-type cooling have shown that these distribu- 
tions retain an approximately log-parabolic sh ape around their 
peaks , with somewhat sharpened curvatures (Massaro et al. 



2006 in particular their Fig. 4). 

The view is widened on using a Fokker-Planck equation 
that includes accelerations and slow cooling of an initial 
mono-energet ic electron injection, as first proposed by Kar- 
dashev]( [l962[ ) and recently computed in detail by Paggi et al. 
(2009a] and the refs. therein). Analytical solutions obtain in 
the form of a power series in the ratio t a / t c of the acceleration 
t a to the cooling time t c defined in Sect. 3.3 The latter reads 
ta/tc = 1.8 X 10" y 2 w/E in cgs units; on using the values 
of w obtained from our fits, this is found to be small out to 
y ~ 10 5 for realistic values of the effective accelerating fields 
E ~ 10 B, in terms of the natural scale B. Then the zeroth 
order solution is the log-parabola, while the first order correc- 
tion steepens the energy distribution at the high end y > 10 5 , 
mimicking a somewhat sharpened curvature that affects the 
SEDs at frequencies well beyond their peaks. 

All that said, it is to be stressed that out to photon energies 
exceeding as 1 GeV emitted by electrons with y > 10 4 , broken 
power laws for the energy distribution (as ado pted, e.g., by 
Donnarumma et aT7||2009[ |Vittorini et al.|[2009] | yield closely 



similar shapes for the SEDs around their peaks; even less sen- 
sitive is the hump of y 2 N(y) also given by the integral in Sect. 



2.1 that contributes the most to the energetics. We conclude 
that the different flare spectra observed in the two blazar sub- 
classes: BL Lacs and FSRQs, intrinsically relate to their dif- 
ferent physical properties during pre-flare conditions, i.e., to 
dearth or plenty of gas surrounding the central SMBH. 

In fact, we have studied in Sect. [3] how this physical di- 
chotomy is reflected in the flaring activity. We find that during 
the flares the BL Lac and the FSRQ follow different patterns, 
indeed diverging ones on the Lj - y p plane. 

We focused on the pattern derived from recent multi-/l ob- 
servations of the two FSRQs, 3C 279 and 3C 454.3. We find 
these flares to run closely vertical on the Lj - y p plane to 
be driven mainly by an increasing flux of the external pho- 
tons, under further growth of the already substantial accretion 
rates. Meanwhile the electron energies hardly increase, be- 
ing constrained by the shortening cooling times in the denser 
radiative field. 

These sources add to three BL Lacs that we previously stud- 
ied (S50716+714, Mrk421 and Mrk501) and complete our 
picture. We interpret their flares in terms of y p increasing un- 
der further acceleration, to yield patterns slanted toward the 
y axis; an effect reinforced by saturation of the power Lj due 
to the limited output extractable from a rota ting SM BH via 
the Blanfor d-Znajek mechanism (Vittorini et al. 2009; Paggi 
|et al.|2009b) . 

So we propose the following pattern for the blazar flares 
analyzed here: flaring spectral changes relate to the pre-flare 
positions of BL Lacs and FSRQs on the Lj - y p plane, specif- 
ically, to their respective branch of the bright blazar strip; this 
is because the constraints that set such positions are retained 
or reinforced during the flares. To wit: bright sources lying in 



the upper branch of the strip flare up almost vertically as ex- 
pected from the cooling constraint to y, while sources in the 
lower right branch move almost horizontally when luminous, 
as expected from the BZ constraint to output given their BH 
masses M # ~ 10 8 H- 10 9 Mq. 

Such a picture may be tested on more sources with Fermi 
and multi-wavelength data; in particular, it will be fruitful to 
study any interlopers between gas-rich powerful FSRQs and 
gas-poor BL Lacs (that is, lower-luminosity FSRQs and the 
LBLs, low-peaked BL Lacs) during their flares on the Lj - 
y p plane in search of any divide or smooth rotation between 
these patterns. This may be the case for the BL Lac sources 
PKS 0537-441, AO 0235+164 and PKS 042 6-380; with their 
weak broad lines and with m ~ 0.1 (see |GMsellini et al.|2009) > 
they may constitute transitional objects between FSRQs and 
BL Lacs. Further investigation of such objects will help un- 
derstanding their nature and their stance in the above picture. 
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Table 1 

Model parameters for 3C 454.3 and 3C 279 states described in the text; R is given in units of 1()' 6 cm, n in units of cnT 3 , B in G, D in units of 10 16 cm, L r and 

Lt in erg s _1 . 



Source 


n 


y P 


r 


R 


r 


B 


D 


fh 


U 


L T 


3C 454.3 (low state) 


5 


421.5 


0.9 


71.9 


16 


0.2 


50.9 


0.4 


1.4 x 10 47 


3.1 x 10 47 


3C 454.3 (high state) 


8 


425.3 


1.0 


61.3 


16 


0.2 


51.2 


1.0 


2.8 x 10 47 


4.6 x 10 47 


3C 279 (low state) 


263 


130.3 


1.1 


4.7 


15 


1.3 


4.4 


0.3 


3.5 x 10 45 


3.4 x 10 46 


3C 279 (high state) 


140 


137.3 


1.1 


6.2 


15 


1.4 


24.3 


2.7 


6.4 x 10 45 


3.9 x 10 46 
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Figure 3. Flaring p atterns of blazars on the Lt - Jp plane, superimposed to the bright bla zar st rip, covered by historical blazar states, with data adapted from 
Celotti & Ghisellini (20081 on taking for y p the peak of their y 2 N(y) distributions (see Sect. |2.l| . Light blue dots represent BL Lacs and light red dots represent 
fSKys. 1 ne arrows snow the patterns of evolution: for gas-poor but luminous BL L acs ( blue arrows) the flare directions turn out to be substantially slanted to 
the right; for FSRQs (red arrows, extended by a factor 3 for better visibility, see Sect. |3.3) we find closely vertical directions. Our actual high and low states are 
reported in detail by Paggi et al. 1 2009b t for BL Lacs, and in Table^for the FSRQs. The black arrow represents the flare pattern of the weaker and more massive 
Mrk501, not affected by constrains on luminosity or electron energy. 



